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■ ABSTRACT 

We present a three dimensional model for the Milky Way fit to the far- 
infrared (FIR) and near- infrared (NIR) data from the COBE/DIRBE instrument 
for galactic latitudes |6| < 30 degrees and to within 20 degrees of the Galactic 
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O^f center. Because of the low optical depth at 240/im, the FIR emission traces the 

distribution of Galactic dust in the Galaxy. We model the dust distribution as 
^ ' due to three components: a warped exponential disk with scale length 0.28 R & 

and a flaring scale height, a spiral arm component with four arms as traced by 
Galactic HII regions, and the local (Orion) arm which produces prominent emis- 
sion features at galactic longitude I ~ 80 and —100 degrees. A Cosmic Infrared 
Background of 1.07MJysr _1 is recovered, consistent with previous determina- 
tions. The dust distribution is then used to calculate absorption in J and K, 
and the stellar emission in these wavebands is modeled with two components: a 
warped exponential disk with a scale length of 0.28 R & and a spiral arm com- 
ponent dominated by two arms. This small scale length is consistent with a 
maximal disk model for our Galaxy, which is inconsistent with the cuspy dark 
matter halos predicted in CDM models. We find different amplitudes for the 
warp in the stars and dust, which starts within the Solar Circle. 

Subject headings: galactic structure, dust, IR, absorption 
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1. Introduction 

The detailed study of the structure of the Milky Way enables astronomers to address 
many of the most important questions in astrophysics, for it is only in the Milky Way that 
we can make detailed studies that enable us to infer the distribution of dark matter, the star 
formation history of the Galaxy, the evolution of the Galaxy's spiral structure and morphol- 
ogy. The COBE satellite's all sky near- and far-infrared maps are a new window into the 
structure of our Galaxy and the beginning of a bridge between optical and radio observations 
of the Milky Way Its potential in revealing the Galaxy lies in the greater transparency of 
the interstellar medium (ISM) to light at these longer wavelengths. The NIR emission arises 
mainly from stars, as at visual wavelengths, but probes greater distances at low galactic 
latitudes, while the FIR emission is an even more powerful probe of Galactic structure as 
the Galaxy is nearly transparent at these wavelengths, much as in the radio regime. The 
FIR emission is from warm dust - the same dust that is responsible for absorption in the 
visual and NIR bands - and therefore uniquely traces an important component of the ISM. 

Due to extinction effects previous studies of Galactic structure based on visual Galactic 
emission, and typically using star count analysis, have been restricted to high galactic lati- 
tudes, and therefore have been relatively insensitive to the non-axisymmetric structure in the 
stellar distribution. However, with NIR data large-scale asymmetry in the Galactic stellar 
distribution begins to be clearly revealed in a Galactic warp and the barred structure of the 
Galactic bulge. The existence of a Galactic warp has been known since the first radio sur- 
veys of the Galaxy, but was not seen in the stellar distribution until Djorgovski and Sosin's 
analysis of IRAS data (Djorgovski and Sosin 1989), and then later by Freudenreich et al. 
(1994) in the DIRBE NIR data. Meanwhile the nonaxisymmetric structure of the Galactic 
bulge was first revealed in the IR by Blitz and Spergel (1991), and later confirmed in the 
NIR from COBE data by Weiland et al. (1994). Since these discoveries both parametric 
(Freudenreich 1996, 1998) and nonparametric methods (Gerhard and Binney 1996; Binney 
et al. 1997) have been applied to COBE data to characterize these nonaxisymmetric features. 

The Galactic bar and warp are now recognized features of our Galaxy, yet another 
endemic nonaxisymmetric feature of disk galaxies is spiral structure. In our Galaxy these 
are perhaps best mapped by HII regions (Georgelin and Georgelin 1976) and CO emission. 
In this contribution we apply parametric models to both FIR and NIR COBE data in order 
to describe the nonaxisymmetric structure in our Galaxy beyond r = O.35i? . 

Other than the studies already mentioned above, the IR and FIR data has also been 
modeled (after the removal of stellar emission in the IR by an assumed stellar emission model) 
by decomposing the emission into two to three components, each associated with a different 
gas phase of the ISM, taken to be correlated to HI, CO, and HII radio fluxes (Sodroski et al. 
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1997, and references therein). Applying a dust emission model to the recovered spectra then 
allows physical state variables of the dust to be inferred. This approach efficiently allows 
the removal of the Galactic foreground emission to recover the Cosmic Infrared Background 
(CIB), and sheds light on dust properties, but it does not yield a detailed map of the dust 
distribution; the decomposition procedure is limited to describing gross radial variations of 
dust properties, in part due to the kinematic distance ambiguity within the Solar Circle. 

In contrast to the above approach we use a parametric model to characterize the spatial 
distribution of the dust, making as few assumptions as necessary regarding microscopic dust 
properties. First efforts in this direction can be found in Spergel et al. (1996) and Davies et 
al. (1997), who both compare the COBE/DIRBE FIR data with the predicted emission from 
axisymmetric models of the dust distribution. Here we employ a more complex model with 
multiple flux density components including not only an (cool) axisymmetric disk, but also 
(warm) nonaxisymmetric features needed to account for major features in the FIR emission 
profile of the Galaxy. We adjust the model parameters by fitting the predicted 240/im sky 
emission to the 240/im Galactic emission as seen in the COBE/DIRBE instrument. 

Once a dust model is constructed using the FIR emission we in turn use it to account for 
absorption in the NIR, and apply a parametric stellar flux density model to the DIRBE J and 
K band data. Like previous Galactic stellar distribution models, we employ an exponential 
disk component, but with the addition of spiral arms. A Galactic warp is applied to all 
components. Our efforts are similar to that of Freudenreich (1998) (hereafter F98), though 
with some important differences, particularly in the derivation of the dust distribution from 
FIR emission rather than from NIR absorption, as in Spergel et al. (1996). We also limit our 
analysis to the Galactic plane (\b\ < 30 deg), and exclude the Galactic center (GC) region 
(|/| < 20 deg), thereby obviating the need to model the Galactic bar and allowing us to focus 
on the nonaxisymmetric structure at galactocentric radii r > .35i? . 

In the following section the data reduction is described while in Section 3 the dust 
and stellar models are detailed. In Section 4 the results of the parameter fitting procedure 
are given, and Section 5 then discusses the uncertainties of the parameters. The final two 
sections are reserved for discussion of our results in light of previous work, and summarizes 
our conclusions. 



2. Data 



In this work, we analyze the 'Zodi-Subtracted Mission Average (ZSMA)' maps produced 
from the DIRBE data by Kelsall et al. (1998) for the 240 pim band in the FIR, and the J 
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and K bands in the NIR. The initial data set consists of a total of 393216 pixels in each 
waveband on a projected cube. Full skymaps of this data are shown in Arendt et al. (1998). 
The data reduction procedures discussed below were done using the UIDL data analysis 
package developed by NASA's Goddard Space Flight Center, as well as additional IDL code 
developed to work in the UIDL environment. 

From the NIR bands the point sources are first removed by applying a median filter to 
the data, without smoothing, to identify and reject point sources (pixels with much higher 
intensities than neighbors) and questionable data points (pixels with much lower intensities 
than neighbors). Using a 5 x 5 pixel window a local median intensity is found at each pixel 
location, and pixels are kept whose intensity is between 0.5 and 1.5 times the local median. 
This procedure was applied to each NIR waveband separately, but any pixel rejected in one 
NIR band was rejected from the other NIR band as well. A local estimate of the uncertainty, 
or error, <t\, is made at each pixel location in each waveband (A = J,K), via the standard 
deviation of the emission of non-rejected pixels in a 7 x 7 pixel window. The estimated error 
for the J band is shown in Figure 1 with respect to galactic latitude. The K band error map 
is similar. 

Unlike the NIR emission, which comes primarily from unresolved point sources, the FIR 
emission is inherently diffuse and possesses structure at all angular scales. As a result the 
emission locally appears as diverse features such as "ridges" and "peaks" . This necessitates 
that a different procedure be used to evaluate the local error. In addition, because the signal 
to noise is much lower in the 240 /zm band than in the NIR bands, pixels are found with 
negative intensities after Zodiacal light removal. These values are obviously unphysical, but 
we are want to simply exclude these data, as doing so would introduce a bias by effectively 
adding a net positive flux. Some smoothing is required, but with care so as to not system- 
atically redistribute flux near the Galactic plane (GP), where the gradient of the emission is 
high. 

At each pixel location a second-order two-dimensional polynomial is fit to the FIR 
emission in a 7 x 7 pixel region. The value of the resulting fit at the central pixel location 
is taken as the smoothed value, while the standard deviation of the residuals in the window 
is used as an estimate of the local uncertainty, <x 2 4o- This procedure introduces minimal 
smoothing with less noise than the data; no negative intensities are present in the smoothed 
data, and the total smoothed flux is equal to the total observed flux to two parts in ten 
thousand. It is to this smoothed data that our model of the the dust emission is fit. Figure 
2 shows a map of <t 2 4o and the ratio of the smoothed to raw data, demonstrating that no 
redistributed flux associated with the GP is present in the smoothed data. 

Here some words concerning the purpose of a local estimation of the error are appropri- 
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ate. 

The local errors are used in the calculation of the merit function used in the parameter 
fitting procedure to be discussed below (see Section 4). With COBE's large beam and its high 
sensitivities, we are in the "confusion" limit where stars just below our detection threshold 
and small scale structures in the ISM are the dominant source of "noise". By using local 
deviations from a local fit or local mean as a measure of the error, an empirical estimate 
is being made of how "noisy" the data is in each direction. This noise is due not only 
to instrumental or measuring error, but also includes those variations arising from (often 
unresolved) structure on a scale smaller than the model intends to describe. In this sense 
"uncertainty" is a much better descriptor than "error". For the NIR bands the primary 
source of uncertainty is unresolved stars, while in the FIR the uncertainty arises from both 
unresolved and resolved structure. In the merit function those pixels with higher error are 
given less weight than those with smaller error, and are thus allowed to deviate further from 
the model. 

Finally, as we are primarily interested here in the distribution of stars and dust in the 
Galactic disk, all data with |6| > 30° are excluded. Also excluded are pixels with longitudes 
within the 20 degrees of the GC, the Orion Nebula, and regions centered on the extragalactic 
sources Andromeda, M33, and the Magellanic Clouds. This leaves a total of 173569 pixels in 
the 240/zm dataset, while after point source removal there are 152371 pixels remaining in each 
NIR band, thus giving a total of 478311 data points for the adjustment of the parameters. 

3. The Model 
3.1. preliminaries 

Of the wavelengths at our disposal, the FIR bands provide the most direct indication 
of the distribution of Galactic dust, whereas the NIR bands reflect the stellar distribution 
and give evidence of the dust distribution through both absorption and emission features. 
We first consider the nature of the FIR measurements before turning to the details of the 
model. 

The 240/im emission, after Zodiacal light has been subtracted, is due to Galactic dust 
and isotropic extragalactic emission. Using the same approach as Spergel et al. (1996), we 
ignore self-absorption or scattering at 240/im as the optical depth is much less than one 
through the Galactic disk. The emission in any particular direction is then proportional to 
the dust column density along the line of sight and the dust's emissivity. Explicitly, the 
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received specific intensity from a line-of-sight due to Galactic FIR emission is 

POO 

I V {T) = / p{s)e v {T) ds, (1) 



where p and e u represent the dust density and temperature dependent specific emissivity 
respectively. However, the actual quantity represented in the skymaps from the DIRBE 
instrument are inferred specific intensities based on the assumption that I u u is constant. If 
the actual spectrum is otherwise then a color correction term is needed to relate the inferred 
intensities to the actual intensities: 

D h = I h K h , (2) 

where is the correction term and ib and D b are the actual and inferred specific intensities 
at the representative frequency b. The correction term 

L 



Kb(T) = I fW h (u) du, (3) 



Wb(u) being the normalized frequency response of the instrument for the b passband (Hauser 
et al. 1999). 

If the intensity arrives from a component with a single temperature, then from equations 
(1) and (2) 

D b = e h (T)K h (T) j p(s)ds = e h (T)K h (T)n, (4) 

with n being the dust column density. If a sum of single temperature components are present 
then 

D h = ^2 thiT^KbiT^rii = ^ h(Ti)ni, (5) 

i i 

where 

^ = e h (T o) K h (T o) (6) 
are relative emissivities, with the (unknown) constant eb(T )Kb(T ) being absorbed into 
the column densities rij. This constant remains unknown as long as T and e„(T) remain 
unspecified. Indeed, if such a model for the emissive medium is assumed (as in the earlier 
work of Spergel et al. (1996)), it is sufficient to define k = 1 for one of the components 
to proceed, letting the absolute temperature and specific emissivity of the medium remain 
unknown. 

In the case where the temperature is spatially varying a more complicated treatment is 
necessary. From equations (1) through (3) 

D h = J W b (i/) du = J p(s) J e u (T)W h (u) duds = J p(s)e h (T)K h (T) ds, (7) 
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where we have defined a new (local) color correction term 




W h {v) du. 



(8) 



This new color correction term is spatially varying inasmuch as the temperature varies spa- 
tially but will be equivalent to the previous color correction term (equation 3) if the tem- 
perature is spatially invariant. 

If it is assumed that the dust emissivity e v oc v a B v (T), where B V (T) is the Planck 
function, then we can define a relative emissivity 



absorbing the constant B h (T )K h (T ) into the normalization of p. 

It is important to point out that there is an inherent ambiguity in scale when confronted 
solely with a 2D intensity skymap, which arises from the integrated flux column density of 
3D emitting structures. For arbitrary distributions of emitting matter, the relative size and 
flux density are mutually indeterminate; a given distribution will result in an intensity map 
identical to another with twice the size and half the flux density. Therefore dimensionless 
units are used in the model by setting R Q = 1, effectively fixing the relative scale of the disk 
and spiral arm components of our model. 

In what follows (r, 0, z) represent galactic cylindrical coordinates, being taken in the 
direction of Galactic rotation and with the Sun lying along = at (R e ,0, Z Q ). The 
transformations used to go from heliocentric galactic coordinates (/, b, s) to galactocentric 
Cartesian coordinates are 



B h (T)K h (T) 
B h (T )K h (T o y 



(9) 



and write 




(10) 



x = R Q — scos(6) cos(/) 
y = s cos(6) sin(Z) 
z = s sin(&) + Z Q . 



(11) 



We are now prepared to describe the model of the dust distribution. 
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3.2. dust emission model 

The inferred specific intensity from Galactic dust emission is modeled as arising from 
three density components, plus an isotropic contribution: 

D 240 d (l, h )= I (Waxi + ^aParm + k + „p\ oc )ds + Q 2 40 , (12) 
JO 

where each component has an associated relative emissivity kj, j = d,a,+,—. Q240 is an 
offset applied to account for isotropic extragalactic flux, i.e. the cosmic infrared background 
(CIB). 

The earlier work of Spergel et al. (1996) described the dust density with an axisymmetric 
exponential disk and implicitly assumed a constant emissivity (i.e. temperature); similarly, 
an exponential component is employed here with the form 

Paxi = p exp(-r//i r )sech 2 (z//i d ), (13) 

where h r and /i d are the scale length and height respectively. However, the scale height of 
the disk is given a linear flair, so that 

7 / \ f h , d + h lid (r - r { ) r > r f 

h d (r) = < 14 
I h , d r < n . 

A hole is effected in the dust disk by taking p axi — > p ax i(?" = 0.5i? Q ) exp(— (r — 0.5) 2 /0.25 2 ) 
for r < O.5i? . No inner dust ring is employed in the model. 

For the disk component we assume a linear radial temperature gradient of — 6.8K/i? Q 
with a central (reference) temperature of T Q = 26K and minimum temperature of 3K. This 
temperature profile approximates that of the dust associated with neutral hydrogen, as found 
by Sodroski et al. (1997). From equation (9) kd(T) is computed, using the fitted polynomial 
of Schlegel et al. (1998) for K h (T), consistent with the assumption that the dust emissivity 
is well described by e v oc v 2 B u (T) (Draine and Lee 1984; Dwek et al. 1997). For the range in 
temperatures assigned to the disk, the emissivity is approximately linear with temperature 
out to 1.5i? (see Figure 3). 

To this exponential disk component two additional single temperature components are 
added, being needed to describe both prominent FIR emission and NIR absorption features, 
namely the spiral arms of the Galaxy and a local feature, described as a spiral arm segment. 
Their respective emissivities are adjusted and relative to kd(T ). 

For the spiral arms the geometry of Georgelin and Georgelin (1976) is adopted, as 
implemented by Taylor and Cordes (1993), who map out four major spiral arms based on 
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the location of HII regions. The density profile across an arm is Gaussian, with a half-width 
w a oc r in the GP, and a flaring scale height h a . For any given position in the Galaxy 
the nearest point of the ith spiral arm can be found numerically, with the corresponding 
minimum distance d^i = 1, ..,4. (The enumeration of the arms follows that of Taylor and 
Cordes.) The contribution of the spiral arm component to the density is taken to be from 
that arm for which the quantity exp — (di/w a ) 2 is a maximum. We then have: 

Parm = Pa#a exp - (di/w a f exp - (z / h a f , (15) 

where p a is the density normalization for the arms and g a describes the radial cutoff of the 
density along the center of the spiral arms, taken to be 

„ („ \ / exp(-(i? a -r m f/rl) R a >r m 

x 1 L a _ ' m j 

R a being the galactocentric radius of the nearest point of the spiral arm. The scale width and 
height of the arms, w a and h a , are functions of R a ; the width is assumed to be proportional 
to galactocentric radius (w a = c a R a ) and the scale height of the arms is given a quadratic 
flair: 

h ( p \ / + ^l,a(^a ~ ?"f,a) 2 #a > r f , a 

W = \ h , a R a < r Xa . (17) 

To the description of the spiral arms we add a final refinement, namely a reduction factor 
/ r on the size (c a and h a ) of the i = 3 (Sag-Car) spiral arm. The need for such a reduction 
factor can be seen immediately from the emission profile in the GP; the inequality of the 
emission peaks at ±50 deg longitude, arising from tangent points approximately equidistant 
from the GC, show that the i — 2 and 3 arms are not equal. 

The inner spiral arms account for emission features within |/| ~ 80 deg near the GP, 
producing peaks in directions corresponding to the tangents of the arms, while the outer 
(Perseus) arm produces a broad emission feature from approximately / = 90 to 180 degrees. 
However, inspection of the FIR emission profiles at low galactic latitudes also shows two 
prominent features that lie at about I = 90 and —100 degrees. These peaks correspond with 
prominent absorption features seen in the NIR data, showing them to be produced by local 
structure. The direction of the features corresponds to the local (Orion) "arm" first seen 
in the distribution of young stars. While responsible for the majority of young stars in the 
vicinity of the Sun, the study of HII regions by Georgelin and Georgelin (1976) showed this 
local "arm" to be a minor feature relative to the major spiral arms in the Galaxy. 

To model this local feature a spiral segment with a gaussian density profile is employed: 

Pioc = p s exp(-d 2 /w s 2 ) (18) 
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where d% — (r — r s ) 2 + (z — Z s ) 2 , with r s = i? s exp(— a s 0) describing the spiral segment in 
the GP with a pitch angle p s = tan(a s ), and Z s being the height of the local arm from the 
warped GP. It was found necessary to place the Sun within a gap, parameterized by the 
azimuths of the gap boundaries, 0i and 02- That the Sun resides in a diffuse region is also 
suggested by other studies of the local interstellar medium (Paresce 1984; Frisch 1995, 1996). 
For computational convenience the local arm is truncated at the heliocentric distance where 
it becomes unresolved within a pixel. Beyond the truncation points and within the gap, d s 
becomes the distance to the nearest end point of the local arm. Finally, the local arm is 
given a different emissivity at positive verses negative longitudes (&+,_), as it clearly shows 
evidence of having different temperatures: in emission it is much more prominent at positive 
longitudes while its NIR absorption features at positive and negative longitudes are very 
similar, and the FIR color ratio (D UQ / D 2 aq) shows it to be significantly hotter at positive 
longitudes. 

To resolve the ambiguity in scale the geometry of this local feature must be fixed with 
respect to the GC. This is effected by fixing one of its geometrical parameters, <f>i, and 
adjusting the ratios (w 8 /(f>i), (Z B /cf>i), and (02/0i)- 111 principle the additional information 
contained in the NIR absorption features associated with this local arm enables this structure 
to be accurately placed. 



3.3. stellar emission model 

We assume that the NIR flux is due to stellar emission, moderated by dust absorption: 

^ + K v p A I v = p*e„, (19) 
as 

k v representing the specific opacity of the dust and e„ the specific emissivity of the stars. 
Integrating over frequency after multiplying by the passband's frequency response function 
W / b( z/ ) ) an d remembering that D h = f I v W^{u)du we arrive at the equation 

dD^ 

+ n h p d D h = p*e h K h , (20) 

where we have used the color correction factor defined in equation (8), though in this 
context it is not a function of temperature. The densities pd and p* are for the dust and 
stars respectively, and the opacity is now a mean opacity: 

fK u I„W h (u) dp 

Kh ~ fi„wMdv ■ (21) 
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If the stellar NIR emission is dominated by a single stellar population, possessing a charac- 
teristic specific emissivity e v (or luminosity function), then eb and K h are spatially invariant. 
We therefore combine all the terms on the RHS of equation (20) into a stellar flux density, 
?7b, and arrive at 

D£° d (l, b) = [ Vh (s) exp(-r b (s)) ds + Q h . (22) 

Jlos 

The mean optical depth r b (s) is found from the dust distribution model and Qb is an isotropic 
offset term, albeit here without astrophysical justification. 

The optical depth r b (s) in equation (22) is defined as 

r b (s) = / n h p A ds . (23) 
Jo 

We then use the following approximation, treating k\> as spatially invariant: 

r b (s) « n h I p d ds' = K V (A h /A v ) I p d ds' . (24) 
Jo Jo 

We adopt the (Aj/Ay) and (A K /Ay) ratios of Rieke and Lebofsky (1985), and leave k v as an 
adjustable parameter. Rather than using the dust column density as given by the model for 
calculating the opacity, we follow Spergel et al. (1996) and use a rescaled density, achieved 
here by applying a line-of-sight scaling factor to only one of the density components of the 
dust model. The scaling factor is determined by requiring that the predicted FIR emission 
using the rescaled density be equal to the observed FIR emission. For a given line-of-sight 
the modeled FIR emission is 

D 24o d = + <224o, where the sum is over the emission 

contributions of the three dust density components. The appropriate scaling factor for 
component j is then 

j. = D 240 - A ~ Q2A0 ^ 

Dj 

For each line-of-sight one of the scale factors fj is chosen and used to rescale that component's 
density, i.e. pj = fjPj, and the total dust density pd used in equation 24 is recomputed with 
this rescaled density. The component chosen to be rescaled is that whose scaling factor results 
in the smallest fractional change in the component's column density, that is, the scale factor 
fj which minimizes |1 — fj\. However, if the modeled intensity from the nonaxisymmetric 
components is at least 10% of the observed intensity then either spiral arms or the local 
arm is rescaled. This bias only has affect in or near the GP where the condition is satisfied. 
It is also required that fj > to be considered as a valid rescaling factor. This rescaling 
procedure assumes that the entire residual in the FIR intensity, D^q — D™ d , is due to the 
modeled dust density deviating from the true density, rather than from deviations in the 
emissivity (i.e. temperature). 
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If more than one stellar population is modeled, then the RHS of equation (20) will consist 
of a sum of components, namely Yli P*,i e b,iKb,i — YliVb,i, where only the normalization of 
each component is waveband dependent. The stellar flux density r?b is modeled using two 
components, an exponential disk and a spiral arm component, both with a sech 2 vertical 
structure. The stellar spiral arms are assumed to have a density proportional to the disk 
component, so that the total stellar flux density is then 

Vh = r& exp(-r/r*) [sech 2 (z/K) + B h g(R a ) exp -(^/<) 2 exp -(z/K) 2 ] . (26) 

As for the dust, di is the distance in the GP to the nearest spiral arm, the half width u>* is 
proportional to i? a , the radius of the nearest point of the spiral curve, h* & is the scale height 
of the arms, assumed to be constant with galactocentric radius, and g(R a ) is a function that 
describes the variation in density along the arm. The parameter (b = J, K) describes the 
relative amplitude of the arms; allowing this parameter to vary with waveband is equivalent 
to assuming that the arms consist of a stellar population that differs from that of the disk 
component. In addition, a stellar cutoff factor is applied for radii greater than the cutoff 
radius r c : 



/cut = exp j , (27) 

applied so that ??b — ► /cut^b for r > r c . 

We explore two different basic geometries for the stellar spiral arms. One is of a log- 
arithmic form, whose phase and pitch angles are adjusted. In this case the n th arm of m 
logarithmic spirals has a radius 

R n = Rq exp(— a0) exp — (n — 1)^ ,n — l,..,m (28) 

where a is the tangent of the pitch angle, p. For these spirals g(R a ) = 1. Both m = 2 and 
m = 4 spiral geometries will be tried in Section 5 against the data. 

The second spiral arm geometry, adopted for our standard model, is a sheared version 
of the spiral model employed for the dust distribution. Under the assumption that the dust 
spirals show the location of star formation fronts in the Galaxy, a young stellar population 
may be expected to drift "downstream" from the arm by an amount described by an offset 
in galactocentric azimuth. If the arms are assumed to have a fixed pattern in the presence 
of a flat rotation curve, then the drift in azimuth is 

^ = l '° T (;-£)• (29) 

where V Q is the circular rotation speed, r the mean age of the population, and Rc the 
corotation radius. The above shear in azimuth is applied to the loci of the points describing 
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the dust spirals. Because of mean color differences in stellar populations of different ages, 
we expect the spiral arms to have different offsets in the different wavebands. Thus we let r 
vary with waveband (henceforth t^), as well as the width of the arms (hence = c^R a ). In 
this formulation the geometry is found by adjusting parameters equivalent to V t^ and Rc- 
For these sheared spirals the function g(R a ) is the same as for the dust, and the reduction 
factor on the Sag-Car arm is applied only to its width. It is this model of the spiral arms 
that we adopt as our standard model, shown in detail in the next section. 

A global warp is added to all components by making the substitution z — > z' in the 
above formulation, where z 1 = z — Z w for the dust components, or z — for the stellar 
components, the function Z w (r, 0) describing the vertical displacement of the warp. For the 
dust 

Z w = K{r-) sin(0 - W ) , (30) 
where W is the phase of the warp, and the amplitude function 

M0 = { ^ w(r " rw)2 (31) 

with r w being the galactocentric distance that the warp starts, and a w an amplitude param- 
eter. For the stars a separate amplitude coefficient a* w is applied, allowing the stellar warp 
to have a different amplitude, similar to F98. 



4. Adjustment of parameters 

In order to determine the set of parameters that yields a best fit to the data we minimize 
a merit function of a Chi-squared form, specifically the sum of the y 2 for each waveband: 



Xb 



E 



/ 7~)mod 7~)obs\ 

{Uh Dh j (32) 



the second summation being over all pixels, and <Tb being the uncertainties estimated in 
section 2. The general nonlinear minimization routine N2FB from the PORT mathematical 
library, developed by AT&T Bell Laboratories, is used to minimize x 2 - 

The determination of the intensity Dj 3 nod (/, b) predicted by a model entails the numer- 
ical evaluation of the line-of-sight integrals in equations (12) for the FIR and (22) for the 
NIR. These numerical integrations are carried out by applying the trapezoidal rule with an 
exponentially increasing quadrature. The exponential convergence of the integrals allow us 
to approximate the infinite range by integrating out to a finite number of effective scale 
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lengths; we integrate to ten effective scale lengths in one hundred steps. Intermediate val- 
ues of the line-of-sight integration through the dust are tabulated to determine the column 
density f Pd(s) ds in equation (24) to calculate the optical depth t(s). To test the numer- 
ical integration, comparisons were made with analytic solutions of an axisymmetric case, 
namely p oc exp(— r/h r — \z\/h z ), at various latitudes along the Galactic meridian (1 = 0, 180 
degrees). The results of these tests show a relative error of 2 x 10~ 4 . 

Unfortunately due to computational limitations and the number of parameters that must 
be employed to describe both the stars and the dust, a simultaneous fit to all wavebands 
was not possible. The fitting of the parameters is thus done in a two step process. In the 
first step parameters describing the dust distribution are adjusted using only the FIR data, 
then the remaining parameters to the NIR data. 



4.1. adjustment to the FIR 

The integration of the dust model along multiple lines-of-sight renders an FIR intensity 
skymap which can be directly compared with the DIRBE 240/im skymap. However, there 
are several ambiguities which the FIR data alone does not allow to be resolved. Already 
mentioned is the ambiguity in scale, which leads us to adopt length units such that R Q = 1. 
A second ambiguity remains in the decomposition of the flux densities into emissivities and 
densities, requiring that one or the other be fixed. For the fit to the FIR emission we specify 
that k- = k a = 1., which are relative to the emissivity of the disk component, already 
defined by it's assumed temperature gradient. (The requirement that the local arm has the 
same density and dimensions at both positive and negative longitudes adds an additional 
constraint that allows k + to be adjusted relative to 

While setting R Q = 1 defines the scale of the model, the Sun's height above the GP, Z Q , 
in principle remains a free parameter to be adjusted. However, there is a near degeneracy 
between Z Q and the other parameters that determine the latitude of the projected emission 
from the local arm. For this reason it was found necessary to fix Z & in the FIR adjustment 
while adjusting the warp parameters and (Z 8 /<f>i). We have set Z & = 0.001875 (15pc/8kpc), 
consistent with other estimates based on NIR data (16 pc, F98; 15.5 pc, Hammersley et 
al. (1995)), though smaller than recent determinations based on optical starcount analysis 
(20.5 ± 3.5 pc, Humphreys and Larsen (1995); 27 ± 3pc, Mendez and van Altena (1998)). 
However, Haywood et al. (1997b) point out that these later estimates are sensitive to the 
adopted scale height of the disk, and should be adjusted downward to ~ 15 pc if a smaller 
effective scale height of the disk is adopted, as they advocate. It's worth noting that if the 
local arm is assumed to be in the GP (i.e. Z s — 0), adjusting Z Q leads to Z Q < 0, contrary 
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to almost all studies of Z & which place the Sun above the GP; this was the primary reason 
for introducing Z s as a parameter. 

A second parameter found to be poorly determined is the inner scale height of the spiral 
arms, h 0jSi . This parameter is constrained by the emission of the inner spiral arm tangents, 
but these make up too small a fraction of the total emission to allow a reliable adjustment 
of /i ,a- We have therefore set /i ,a = O.Oli? . 

A total of twenty two parameters are adjusted in the FIR fit: There are four adjusted 
parameters associated with the disk component (p , /i ,d, ^i,d? K), seven associated with the 
spiral arms (p a , c a , r fja , r m , r a , f T ), seven associated with the local arm (p s , R s , p s , (w s /(f)i), 
(Z s /(pi), (02/0i), k + ), three associated with the warp (r w ,a w ,0 w ) and the offset Six 
parameters remain fixed: 0i, rf, /io, a , Z Q and the emissivities k- and k a . The resulting pa- 
rameters of the preliminary dust fit are given in Table 1, expressed for convenience in units 
where R Q = 8 kpc . 

We leave the general discussion of the implications of the parameters for later (Section 
6), but here present the resulting skymap of the FIR emission produced by the model as 
compared with the observations. Figure 4 shows the observed and modeled 240/im sky 
emission, and a map of their relative difference, (D ohs — jj^^yi) ^^ Such grey scale maps 
give an overall impression of the emission, but do not show the variations in the GP where 
most of the details of interest are located. To show a more direct comparison between the 
modeled and observed signal, emission profiles in the GP and other galactic latitudes are 
shown in Figures 5 and 6. In these emission profiles it appears that the model has higher 
residuals at higher latitudes. However, this is an effect of the logarithmic scale used; Figure 
7 shows that the residuals are in fact smaller at higher latitudes. 

The emission profiles show that the major emission features are reproduced. In the GP 
the emission peaks at the tangents of the spiral arms (/ ~ ±30, ±50 and — 80deg), and the 
peaks from the local arm (at I « 80 and — lOOdeg) are correctly placed. In the skymaps 
the local arm shows itself as two bright spots with considerable extent in latitude in nearly 
opposite directions of the sky. The Galactic warp is evident in the skymaps as an asymmetry 
in the Galactic emission at positive versus negative latitudes. An emission feature due to 
the outer Perseus arm is evident as a broad feature from approximately / = 80 to — 140 deg 
in Figure 6 for latitudes b > 5 deg. This arm is seen primarily at positive latitudes due to 
the Galactic warp, and its large extent in latitude as compared to the inner arms is what 
necessitates a flaring in the scale height of the spiral arms. The warp also accounts for the 
small negative deviation in latitude at I ~ —80 deg, a feature associated with dust within 
the Solar Circle. 
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4.2. adjustment to the NIR 

A final adjustment to the NIR data yields the remaining parameters of the stellar 
distribution and resolves the remaining uncertainties of the dust distribution; the additional 
information provided by NIR absorption, proportional to the dust column density, allows 
the decomposition of the dust densities and FIR emissivities. For example, fixing all other 
parameters of the spiral arms, the FIR flux column density from the arms arriving from 
any particular line-of-sight is simply proportional to the product (k a p a ). The decomposition 
is effected by adjusting p a and varying the emissivity k a so as to keep the product (A; a Pa) 
constant, thus preserving the modeled FIR signal. Similarly for the local arm, the column 
flux density is fixed while p s is adjusted. 

In the fit to the NIR the parameter Z Q is adjusted, while the distance to the warped 
plane, Z & — Z w (x Q ), is preserved by appropriately varying W , leaving the parameters a w 
and R w fixed to the values found in the FIR fit. Meanwhile, the amplitude factor for the 
stars, a^, is adjusted to allow the stellar warp to have a different amplitude than that found 
in the dust. 

In summary there are nineteen adjusted parameters in the adjustment to the NIR data: 
five parameters describing with the disk component (rjj K , r*, h*, r c ), seven parameters as- 
sociated with the spiral arms (-Bj 5 K; Cj K , V tj 5 k, -Rc)) three parameters associated with the 
dust and it's absorption (Kv,p a ,p s ), the stellar warp amplitude parameter, a^, the vertical 
coordinate of the Sun, Z Q , and the offsets <5j,k- The stellar model possesses only a single 
parameter that is explicitly fixed, the scale height of the stellar spiral arms, set to be equiv- 
alent to that of the dust (/i* = h 0ja ). However, several assumptions are made with regards 
to common geometry between the dust and the stellar distributions, such as the radius that 
the warp starts. 

The parameters resulting from the NIR adjustment are given in Table 2. The NIR 
skymaps of the DIRBE instrument and of the model are given in Figures 8 and 9, as well as 
the relative differences between the data and the model, (£)° bs _£) m ° d ) / /£)° bs . The skymaps of 
the observations and model look very similar, the deviations only being obvious in the relative 
difference maps. These show that the largest deviations are in the GP within 30 deg of the 
GC. Again, to show in more detail the concordance and deviations between the observed and 
the modeled emission, especially in and near the GP, emission profiles at various latitudes 
are shown in Figures 10, 11 and 12. 

The observed NIR emission profiles are much smoother than those seen in the FIR, and 
in contrast to the FIR the axisymmetric component dominates the emission profiles at all 
latitudes. The profiles deviate most from an ideal axisymmetric profile as one approaches 
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the GP, especially in J, and this is largely due to absorption. The spiral arms are much less 
evident than in the FIR, and only provisionally identifiable in the K band (Drimmel 2000). 
Indeed, differentiating between emission and absorption features in the NIR is problematic, 
and it has been suggested that deviations from axisymmetry in the NIR can be attributed to 
absorption effects alone (Kent et al. 1991). The two absorption features most clearly evident 
are due to the local arm, seen at / ~ 80 and — lOOdeg. 

In the modeled emission profiles the effect of using a rescaled dust model to calculate 
extinction is immediately obvious, introducing fluctuations even on small scales. Incidences 
can be found where the rescaling introduces spurious features, but more often it allows the 
model to achieve a correspondence with the data that could not be obtained otherwise. This 
is particularly evident at low latitudes in J, and shows the importance of absorption at these 
latitudes in shaping the profiles. Nevertheless, emission from the spiral arms is important 
in the GP, though we will have to make comparisons with axisymmetric models in Section 
5 to make this demonstrably obvious. 

One nonaxisymmetric feature that is evident in the NIR is the Galactic warp, though 
less obvious than in the FIR. It reveals itself as an asymmetry in the emission profiles, seen 
most clearly in Figure 12 at |6| = 5deg; the profiles are skewed toward positive longitudes 
at b = —5 deg, and toward negative longitudes at b = 5 deg. This requires a local tilt in the 
stellar distribution with respect to the 6 = plane, accounted for here by the warp starting 
within the Solar Circle. 

The relative difference map in K suggests the presence of structure not accounted for 
by the model. Aside from the deviations in the GP mentioned above, broad bright patches 
in the difference map in K can be seen at positive and negative longitudes. That at positive 
longitudes, roughly located at I — 90 deg, is distinctly above the GP, and may be stellar 
emission associated with the local (Orion) arm. The bright patch at negative longitudes 
is less obvious, centered on the GP and is found within I = —90 deg, thus not correctly 
positioned to be associated with the local arm. Additional features seen in the difference 
map in J are associated with over or under estimated absorbtion. 

4.3. General features of the model 

To concisely summarize the general features of the standard model a "bird's-eye-view" 
of the Milky Way is presented as a surface density map for the dust in Figure 13, and as a K 
band surface brightness map in Figure 14. The center of the maps are determined from an 
extrapolation of the model, which is itself specifically constructed only for r > .34i? , thus 
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the picture presented for r < .34R Q is necessarily incomplete. In particular, no structure 
associated with the Galactic bulge is represented. Also, the spiral arms are incomplete on 
the far-side of the Galaxy due to the adopted spiral geometry based on the observed HII 
regions. The relative strength of the spiral arms in the dust and the stars is shown in Figure 
16 as the arm-interarm (flux) density contrast in the GP. However, the difference in the 
scale heights of the disk and spiral arms mitigate the surface brightness/density contrast 
between the arm and interarm regions; the arm-interarm ratio in the NIR surface brightness 
is 1.2 and 1.32 for J and K respectively. Figure 15 shows the scale heights of the various 
components. 

In addition to the major spiral arms the dust surface density map shows the smaller 
local arm in the region of the Sun. Though the local arm produces prominent FIR emission 
features in our sky, from an extra-Galactic perspective it is revealed to be a minor feature, 
as pointed out by Georgelin and Georgelin (1976). Indeed, "arm" is perhaps a misnomer. 
Similar structures are undoubtedly found throughout the Galaxy in the form of spurs and 
bridges between the main spiral arms; the local arm is prominent and resolved in our FIR 
sky only by virtue of its vicinity. 

Our flux density model implies an extinction corrected K magnitude of —23.79 for the 
Milky Way, assuming R Q = 8. kpc, of which the spiral arm component contributes 7 percent 
of the total luminosity. However, our spiral arm model is incomplete for the side of the 
Galaxy opposite the Sun due to the lack of observed HII regions. Assuming that the spiral 
arm model is complete over 3/4 of the Galaxy, the actual contribution of the spiral arms 
to the total flux would be 10 percent, resulting in a K magnitude of —23.82, which does 
not include light from the Bulge. The bulge/disk luminosity ratio has been approximately 
estimated as 0.2 (Gerhard 2000), which results in a final K magnitude of —24.02 for the 
Milky Way, consistent with the earlier estimate of Malhotra et al. (1996). 

5. Uncertainties and alternative models 

In this section we explore the sensitivity of the adjusted parameters on the choice of 
fixed parameters, thereby estimating systematic errors, and test the relative importance of 
specific features of the model that are not obviously essential. In particular various spiral 
models for the stars are considered in an attempt to reproduce the observed NIR emission 
profiles in the GP. 
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5.1. uncertainties 

Four of the fixed parameters in the adjustment to the FIR, rf, /io, a , Z & ), as well 
as the assumed temperature gradient of the disk component, are each varied in turn from 
their standard values, and the FIR adjustment redone. The resulting values of the adjusted 
parameters are shown in Table 3. In all these fits R s , the radius of the local arm, was adjusted 
to within 0.1 percent of 1-R©, and has thus been left out of the table. The resulting x 2 of 
these fits fell within 1 percent of that of the adopted standard model, with the exception 
of the 0i = 3deg and Z Q = 10 pc models whose x 2 deviated by less than 3 percent. The 
variance of the resulting parameters are used as estimates of the uncertainties given in Table 
1. 

From Table 3 notable correlations between Z Q and adjusted parameters are worth men- 
tioning. The strongest is that between Z & and the linear flare parameter hi of the dust disk, 
and we have excluded these models in the calculation of the uncertainty of hi. There are also 
strong correlations with Q240 an d the spiral parameters c a and hi^. For the local arm many 
of the parameters are correlated to the value of 0i. In the case of the density and emissivity 
this is expected, however the initial intention of the "geometrical" parameters (p2/4>i, w s /<f)i, 
and Z s /(pi was to allow an eventual adjustment of 0i alone in the adjustment to the NIR 
data to arrive at the relative placement of the local arm. If this were possible then these 
geometrical parameters would be independent of the choice of the initial 4>i- However, this 
is not realized. For this reason 0i is kept fixed in the NIR adjustment. 

To estimate the uncertainties in the stellar parameters the ten parameter sets derived 
from the FIR adjustments shown in Table 3 were in turn used for a suite of NIR adjustments, 
allowing an estimate of the uncertainties in the stellar parameters to be made. In contrast 
to the dust model, the stellar model only has one fixed parameter, the scale height of the 
stellar spiral arms, and this parameter is effectively varied in this suite of adjustments because 
hi = /io,a- The resulting sets of parameters are shown in Table 4. 

Inspection of the disk parameters for the different models shows evidence of clustering 
about two solutions, probably a result of two minima in x 2 space, one with an average 
scale length of 2.56 kpc (3 models), and the other with an average scale length 2.19 kpc (for 
R & = 8 kpc). The standard model belongs to the latter group and the uncertainties for 
the disk parameters given in Table 2 is from the variance of the seven models belonging to 
this shorter radial scale length solution. The x 2 only slightly distinguishes between these 
two solutions, the longer scale length solutions having an average x 2 2 percent lower than 
that of the standard model. The variance from all ten models is given in parenthesis. The 
other parameters do not seem to reflect this clustering, with the exception of n v . Excluding 
the three longer scale length models reduces the uncertainty of this parameter to 0.0025. 
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Meanwhile the uncertainties for the dust densities and emissivities are found after excluding 
the 0i = 3 deg and Z Q = 10 pc models, which have spiral arm emissivities that are more 
than three standard deviations from the mean. Again, the variances from all ten models is 
given in parenthesis. 

The magnitude of the relative uncertainties are in most cases between 10 and 15 percent, 
though are higher for the parameters associated with the spiral arms. The data in or near 
the GP, where most of the information resides for the spiral arm parameters, makes up a 
relatively small fraction of the x 2 ? with the result that these parameters are not as well 
determined. Possibly another choice of merit function more sensitive to the GP emission 
would render a better determination of these parameters. 

5.2. Alternative models 

Besides the above models we have also run a suite of models to test the sensitivity of our 
results with respect to the data selection and specific features in the models. Adjustments 
to data restricted to smaller latitudes < 20 deg), or with a model lacking a central 
hole in the dust disk, or with a radial cutoff in the dust disk imposed at 1.75i?©, each 
produced estimated parameters that coincided with the standard model well within the given 
uncertainties. However, the adjustment is sensitive to varying the cut in galactic longitude, 
as including emission within 20 deg of the GC will include emission from structures not 
described by the model, such as the Galactic bulge, and excluding emission within 30 deg of 
the GC would exclude important arm tangents needed to contrain the spiral arms. 

Other alternative stellar models worthy of discussion are given in Table 5. A model 
with a flair in the stellar disk, starting at the same radius as the flair in the dust disk, was 
fit to the data. The resulting flair parameter was quite modest, 6.6 pc/kpc, less than half 
that found for the dust disk, but leads to a moderately larger scale height. No significant 
improvement of the x 2 is seen (< thus no positive evidence for such a structure can be 
inferred. Removing the disk cutoff results in a significant change of only a few parameters, 
while the x 2 was only slightly larger than that of the standard model; this is another feature 
of the stellar disk which is not well constrained. However, the necessity of the NIR offsets is 
confirmed. In this case x 2 increases more than 17 percent when the offsets are set to zero, 
but most of the parameters are not significantly different than the normative values. 

A model in which the dust emissivity of the dust spirals was fixed to unity was adjusted; 
though this has a significant effect on the amplitude of the stellar spirals, other parameters 
are not affected, nor is the x 2 ver Y different from the standard model. This insensitivity 
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shows that the dust density (and emissivity) associated with the spiral arms are not well 
constrained, but likewise neither are the other parameters affected by this uncertainty. An- 
other adjustment was done with the scale height of the stellar spirals set to twice that of the 
dust spirals. This modification significantly changes some of the parameters, including the 
disk scale length, but is not favored as the total \ 2 is 3 percent higher than for the standard 
value. This and other experiments confirm that a small scale height is favored for the stellar 
spirals. 

In order to discern the nonaxisymmetric structure in the Galactic disk from the NIR 
emission, four different models were adjusted against the data, three spiral models and an 
purely axisymmetric model. Figures 17 and 18 show the GP emission profiles for longitudes 
\l\ < 90 deg of the four models, and Table 6 gives their parameters and x 2 in the GP. The two 
logarithmic spiral models use a single amplitude parameter B for both wavebands, consistent 
with the hypothesis that the spiral structure has the same stellar population as the disk. 
Meanwhile the sheared spiral model is formulated under the assumption that they primarily 
consist of young stars born from the spiral arms traced by the dust, and are thus given 
parameters that differ from one waveband to the other. As can be seen in Figures 17 and 18, 
the m = 2 log spiral model and the sheared spiral model perform about equally well, while 
the nonaxisymmetric and m = 4 models fail to produce emission apparently associated with 
the spiral arms, particularly in the directions to the tangents of the Scutum arm, I ~ 30 and 
-50 deg. 

At positive longitudes no clear evidence of the Sag-Car arm is present; the axisymmetric 
and m = 2 log spiral models reproduce the emission profile at Z > 40 deg as well as the other 
two spiral models, though the former two models do not possess the Sag-Car spiral arm. 
However, a closer comparison of the m = 2 log and sheared spiral models suggest that 
the Sag-Car arm is indeed present, at least at negative galactic longitudes; the parameters 
common to the m = 2 log and sheared spiral model are similar, though the very large pitch 
angle of the m = 2 log spiral produces a geometry significantly different from that of the 
sheared spiral model (see Figure 19). This high estimate of the pitch angle is most likely 
due to an effort in the model to account for emission at longitudes / < —50 deg, which is 
more prominent in J. Indeed, adjustment of an m = 2 log spiral model to the K band alone 
gives a smaller estimate of the pitch angle of 15.6 deg (Drimmel and Spergel 2001). All four 
of the models are inadequate in reproducing emission for \l\ < 30 deg. 

In fitting the four above models to the NIR data it was necessary to fix the emissivity 
of the spiral arms in the dust for the axisymmetric case. Decoupling of the emissivity and 
density of the dust spiral component was successful for the other three models, though 
questionably reliable for the m = 4 model that gave an emissivity much higher than the disk 
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component. 

6. Discussion 

6.1. axisymmetric structure 

The scale length of the axisymmetric component of the dust distribution is here found 
to be 0.28-R Q , which is shorter than that of F98 (O.37i? ), and significantly shorter than that 
of either Spergel et al. (1996) (0.48i? Q ) or Davies et al. (1997) (O.62i? ). Apparently a more 
extended dust component is not needed when the distribution of the warmer (nonaxisym- 
metric) components are taken into account. Other models of the Galactic dust distribution 
are based on correlating the dust density with the gas (hydrogen) surface density (e.g. Ortiz 
and Lepine (1993)). Such models point out that an exponential disk model is not appropriate 
in the central regions of the Galaxy (r < O.5i2 ), where the gas surface density shows a hole. 
Nevertheless, a characteristic scale length can be assigned to the outer regions, and these 
show a much wider range than those stated above due to uncertainties in the metallicity 
gradient and the CO:H 2 ratio. The latter is usually taken as constant with galactocentric 
radius, but more recently it has been argued that this ratio may vary considerably with 
galactocentric radius (Sodroski et al. 1997). 

A feature seen in the HI distribution that we have included in our model is a flair in 
the disk scale height. Malhotra (1995), using a Gaussian vertical density profile for the HI, 
finds a scale height that flares from ~ 100 pc to ~ 220 pc going from 0.5 to 1.0 R Q . (Here 
and in what follows, quantities expressed in parsecs are derived on the assumption that 
R Q = 8kpc.) We find a smaller initial scale height and more modest gradient, the scale 
height increasing from 134 pc to 188 pc at the Solar Circle. These values are comparable to 
the constant scale height of F98 (152 pc), but smaller than Davies et al. (1997) much higher 
value of 470 pc. A possible source of bias in our determination of the dust scale height is 
the assumption of no vertical temperature gradient, which if present would result in the 
scale height of the dust being under estimated. However, we expect the vertical temperature 
gradient to be small for the dust disk component, as the scale height of the stars heating the 
dust is significantly larger than that of the dust. 

The isotropic offset parameter Q240 corresponds to the CIB. Our value for the CIB, 1.07± 
0.15 MJy sr _1 , is in agreement with other determinations (1.17 ± 0.53 MJysr -1 , Finkbeiner 
et al. (2000); 0.91 ± 0.15 MJy sr _1 , Lagache et al. (1999); 1.09 ± 0.20 MJysr" 1 , Hauser et al. 
(1998)). All estimates agree when an additional systematic error of ±0.20 MJy sr _1 is taken 
into account. 
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Our estimate of the stellar radial scale length of the Galactic disk (O.28i? ) is shorter 
than older determinations from NIR data (0.38 R Q , Kent et al. (1991)) . We discuss the 
implications of this result in the conclusions. Our determination of the stellar and dust 
scale lengths are most sensitive to their distribution between .5 and 1 R & . Our model is 
not sensitive to a cutoff in the dust density, though we do find a cutoff in the stellar light 
distribution at approximately 10.5 kpc, which is comparable though smaller than that found 
by (Ruphy et al. 1996) from a starcount analysis of DENIS data. However, it is still larger 
than would be expected from the scale length, as compared to external disk galaxies (Pohlen 
et al. 2000). 

Traditional estimates of scale heights for Galactic stellar populations varies from ap- 
proximately 90 to 390 pc, and is commonly correlated with absolute magnitudes, stellar 
type or age in stellar distribution models built to reproduce starcount data (Bahcall and 
Soneira 1980; Ortiz and Lepine 1993; Haywood et al. 1997a). Meanwhile luminosity den- 
sity models which reproduce NIR emission on the sky, like our own presented here, have 
reported scale heights of 247 pc (Kent et al. 1991), 276 pc (Spergel et al. 1996) and 334 pc 
(F98). Our estimate of 282 pc for the scale height falls amoung these values. However, scale 
height estimates do not lend themselves to direct comparison as scale lengths do, because 
the formulation of the vertical density profile does not enjoy a universal consensus as does 
the radial profile. In the past vertical exponential profiles were most common, even after 
it was shown that a sech 2 variation describes an isothermal population in the infinite plane 
approximation. Recently this situation has changed with more sophisticated models that 
employ dynamical constraints that have called into question the traditional exponential ver- 
tical profiles (Haywood et al. 1997a,b). Such models do not possess a simple vertical profile, 
but only approach an exponential far from the GP, similar to a sech 2 profile. However, this 
common behavior among the various proposed profiles does not assist us in making compar- 
isons via asymptotic scale heights, as the integrated luminosity is mainly determined by the 
density variation near the GP. 

We note at this point in the discussion two possible sources of biases in our determination 
of scale heights. The first is the assumption of a constant dust temperature with respect to 
z and it's possible affect on the estimation of the dust scale height, mentioned above, which 
if present would cause an underestimation of the stellar scale height. The other source of 
potential bias is the removal of luminosity from point source removal, which preferentially 
takes place at higher galactic latitudes where the median background light from unresolved 
stars is lower. This effect could lead to an underestimate of the vertical scale height. The 
isotropic term in the emission model may correct for this to some degree, but in the future 
it will be desirable to make a more accurate correction for this effect. 



6.2. nonaxisymmetric structure 



Similar to the results of Freudenreich (Freudenreich 1996, 1998) we find a warp that 
starts within the Solar Circle with larger amplitudes in the dust than in the stars. At 
r = lOkpc our stellar warp has an amplitude of 0.25 kpc (R Q = 8kpc), while F98 finds 
0.18 kpc (his Model S). From FIR emission we find a dust warp 2.7 times larger, whereas 
F98 finds an amplitude that is 1.8 times larger as inferred from NIR dust absorption; his 
dust warp amplitude agrees well with the estimated amplitude of 0.3 to 0.4 kpc that has been 
given for the HI (Burton and Hartmann 1994). Our warp amplitude in the dust does agree, 
however, with the warp inferred from the OB stellar distribution (Smart et al. 1998). In the 
NIR the most important effect of a warp that starts within the Solar Circle is to induce a 
local tilt (tan# = /i w (i? Q )/i? ), as evidenced in the NIR emission profiles at low galactic 
latitudes (Section 4.2). For the stars we find a local tilt of 9 = 0.2 deg while Freudenreich's 
warp renders a local tilt of 9 = 0.5 deg. A tilt of the stellar distribution with respect to the 
conventionally defined b = plane has also been noted in NIR data by Hammersley et al. 
(1995), though they propose that a global tilt of the entire stellar disk with respect to the 
6 = plane is responsible for this feature. The motivation for this alternative model is that 
radio observations seem to indicate the Galactic warp starts beyond the Solar Circle, though 
this is not well contrained by the radio data. Evidence for a warp being present in the stellar 
disk is also found in local stellar kinematics, as seen in Hipparcos data, though it is more 
consistent with a warp starting at or beyond the Solar Circle (Dehnen 1998). While the 
warp both here and in F98 starts within the Solar Circle, this is not a commonly accepted 
feature of the warp. In any case our warp starts at a significantly larger galactocentric radius 
(~ O.85i? ) than that of F98 (between 0.5 and O.56i? ), though this difference may only be 
a consequence of the different functions used to describe the warp amplitude. 

With regards to the spiral arms, it was found that an adopted map of HII regions was 
sufficient to describe the location of peaks in the FIR emission features associated with the 
spiral arm tangents, though it was necessary to introduce a reduction factor on the Sag-Car 
arm. This assumed geometry for the spiral arms is consistent with a four arm model of the 
Galaxy, with a pitch angle of approximately 12.5 deg, consistent with other spiral tracers, 
radio and pulsar data (Vallee 1995). We point out that the dust density associated with 
the arms is not well constrained, as strong absorbtion features from the spiral arms are not 
present in the NIR to assist in the decomposition of the FIR flux density. Furthermore, 
there may be significant temperature gradients associated with the arms as the dust here 
may be primarily heated by young OB stars. By using a single emissivity the widths and 
scaleheights of this component may be significantly underestimated. However, for the same 
reason that the dust density in the arms are not well determined, this uncertainty does not 
strongly affect the estimation of the stellar distribution parameters. 
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Evidence for the spiral structure in the NIR is less evident due to extinction, but im- 
portant in and near the GP. Three different geometries were attempted, the geometrical 
parameters of each model being constrained with the NIR data. Our standard model is 
a sheared version of the same four arm spiral model used for the dust. Spiral models con- 
structed to describe the distribution of NIR point sources have also utilized four arm models, 
though of a logarithmic form (Wainscoat et al. 1992; Ortiz and Lepine 1993). One impor- 
tant difference between the spiral model presented here and these previous models is the 
reduction factor on the size of the Sag-Car arm. The necessity of such a reduction factor 
is demonstrated by the failure of the m = 4 logarithmic model in which all four arms are 
treated equally, and it's effect is consistent with the conclusion that two arms dominate in 
the NIR. This interpretation is further supported by the comparable success of a purely two 
arm spiral model. 

We find that the spirals are stronger in K than in J, which also agrees with observations 
of other spiral galaxies in these wavebands (Grauer and Rieke 1998). The amplitude of our 
spirals, however, are smaller than those seen in most spiral galaxies; Rix and Zaritsky (1995) 
report fractional azimuthal variations in the K surface brightness, (S max — S min )/S min , to be 
of the order unity for their sample of spiral galaxies, while we find here that this quantity 
is equal to (^fnh* & /2h*)B K = .32. Also, the fraction of total light in K from the spirals, 
approximately 10%, is smaller than that seen in most other galaxies (Seigar and James 1998). 
Our weaker spiral arms may be due to their scale height being underestimated relative to 
that of the disk. Such an underestimate could be due to K supergiants dominating the spiral 
flux density in the GP, or be a consequence of underestimating the dust column density 
associated with the arms. Our spiral model also renders an estimate of the corotation radius 
of the Milky Way, which is smaller than other estimates, but which gives a ratio Rc/^* = 2.9 
consistent with determinations of other spiral galaxies (Grosbol and Patsis 1998). 

It is interesting to note that our spiral model is consistent with the dynamical model 
of the Milky Way's spiral arms constructed by Amaral and Lepine (1997), after a trivial 
90 deg rotation, which has four arms but with two arms dominating. Additional evidence of 
spiral arms comes from recent results on the star formation history of the solar neighborhood 
(Rocha-Pino et al. 2000; Hernandez et al. 2000), that indicate an intermittent or episodic star 
formation rate for the Galactic disk. A periodicity of ~ .5 Gyr is suggested by Hernandez 
et al. (2000), and they point out that this periodicity is consistent with crossing a two 
arm spiral. Assuming £l a = 25kms~ 1 kpc -1 , our estimate for the corotation radius, Rc = 
0.83 ± O.O5i? gives a range of star formation periodicity from 0.5 to 0.9 Gyr for a two arm 
geometry. 



-26- 



7. Summary 

We fit joint models for the Galactic dust and stellar distributions to the COBE FIR and 
NIR emission. Our dust model has 22 parameters adjusted to fit 173,569 data points, while 
our parameteric model of the stellar emission has 19 adjusted parameters that fits the large 
scale NIR emission features in 304,742 data points. 

Our model of the Milky Way has several intriguing results: 

(1) We find a small scale length for the stars in the Milky Way disk. Our estimate of the 
stellar radial scale length of the Galactic disk (0.28i? Q ) is shorter than older determinations 
from NIR data (0.38 R & , Kent et al. (1991)) . It is, however, consistent with recent studies 
using NIR data, such as IRAS point sources (0.33 R & , Ortiz and Lepine (1993)), DENIS 
data (0.27 R & , Ruphy et al. (1996)), the Two Micron Galactic Survey (0.25 R & , Porcel et 
al. (1998)), and earlier determinations based on the DIRBE data (0.35 R Q , Spergel et al. 
(1996); 0.31 _R , F98). It also agrees with analysis of local stellar kinematics (0.29 - 0.36 R Q , 
Fux and Martinet (1994); 0.2 - 0.34 R Q , Bienayme and Sechaud (1997)) including recent 
use of Hipparcos data (0.29 - 0.33 R Q , Dehnen and Binney (1998); however 0.19 - 0.25 -R© 
in Bienayme (1999)). Recent studies at visual wavelengths are beginning to converge on a 
shorter scale length as well (0.21 - 0.36 R Q , Ojha et al. (1996), 0.36 R Q , Gould et al. (1997)), 
though it is yet far from unanimous (0.47 - 0.94 R Q , Mendez and van Altena (1998); 0.43 - 
0.56 R & , Ng et al. (1995)). We note that De Jong (1996) finds that the NIR scale length is 
20 percent smaller than the optical scale length in external galaxies. He argues that this is 
the signature of inside-out galaxy formation: the outer regions of spiral galaxies are younger 
than the inner regions. This small scale length has a number of important implications for 
the Galactic mass distribution: (i) it implies that maximal disk models are a good fit to the 
Galactic rotation curve (Dehnen and Binney 1998) and implies a low central density for the 
dark matter halo which contradicts CDM simulations(Sellwood 2000); (ii) can increase the 
microlensing optical depth towards the Galactic bulge and will produce more long duration 
events (Binney et al. 2000; Sackett 1997); (in) it implies that the Galaxy is significantly 
smaller than our neighbor M31 (Hiromoto et al. 1983). 

(2) We find that a two arm spiral structure dominates the NIR nonaxisymmetric emis- 
sion, supporting the earlier simple analysis of Drimmel (2000). What stellar population 
dominates our spiral arms in the near-IR? Diffuse NIR emission associated with the spiral 
arms will have contributions from both young stellar populations, such as K supergiants, and 
from the old disk population if a spiral density wave is present. The geometry from these 
two types of populations will not necessarily be the same, as new stars are born from gas 
subject to hydrodynamic forces. Our standard model of sheared HII arms, that are redder 
than the disk, is consistent with the assumption of the spiral emission being produced by 
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young K-supergiants, while a m = 2 logarithmic spiral model with a more open structure 
and the same color as the disk is consistent with a density wave. The success of these two 
geometries may indicate a need to include both types of spiral emission. The small scale 
height of the arms favors the hypothesis of a young population dominating the spiral emis- 
sion, though the vertical profile of a density wave perturbation in the stars is not known. In 
external galaxies with active star formation, young stars dominate the spiral arms even in 
K band, while in other galaxies the K band light is tracing the older stars and the stellar 
mass (Rhoads 1998). Once the 2MASS data is available for our spiral arms, we will be able 
to address this question in our Galaxy. 

(3) The Galactic warp is here found to start within the Solar Circle and to have different 
amplitudes for the dust and the stars. If this is indeed the case it is an important clue as to the 
nature of the warping mechanism, suggesting that hydrodynamic or magnetohydrodynamic 
forces are important. It may also suggest that the Galactic warp is not a long-lived feature. 
One mechanism that would cause a short-lived warp is an interaction with one or more 
companions; this could cause a different response in the gas and stars of the galaxy. Both the 
Magellanic Clouds (Weinberg 1995, 2000) and the Sagittarius Dwarf (Ibata and Razoumov 
1998) have been suggested in this context. 

However, there is an important caveat to this last result, namely that the only large- 
scale vertical distortion present in the Galactic disk is in the form of a warp. That this may 
be inadequate is suggested by several points. First, the amplitude of the warp in the dust is 
apparently inconsistent with radio data of HI. Secondly, the phase of the warp is found here 
to be less than 1 deg, a result also found by F98. Assuming that the location with respect to 
the Galactic warp is not a determining factor for the existence of observers, the probability 
of our being this close to the line-of-nodes (20 w /7r) is less than 1 in 100. The minimizing 
of 4> w effectively maximizes the local tilt for a fixed warp amplitude; the presence of other 
distortions could produce additional tilt in the dust than what would be produced by a warp 
alone. In addition there are several features to the adjustment of related parameters which 
may be pointing to an inadequacy, namely the inability to adjust Z Q in the initial FIR fit 
and the need for a displacement specific to the local dust feature. 

These points suggest the alternative interpretation that, in addition to a global warp, 
there are small-amplitude oscillations affecting the local structure of the gas and dust, that 
are not described in the model. Evidence for vertical displacements within the Solar Cir- 
cle can be found in the COBE data, particularly with regards to residuals at the spiral 
arm tangents. Hammersley et al. (1999) has also noted evidence that the spiral arms show 
displacements out of the GP in NIR point source data. It has been noted in numerical exper- 
iments that small amplitude vertical displacements in the disk could result from oscillations 
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excited by Galactic satellites (Edelsohn and Elmegreen 1997), though displacements specific 
to the spiral arms would have to be of a hydrodynamic nature. 

We are entering a golden age for Galactic astronomy. Near-infrared and far-infrared 
observations are revealing structure that is hidden in the optical. Our analysis here shows 
that there is a wealth of Galactic structure information in COBE two dimensional data. This 
data will soon be complemented by a 2MASS inventory of the bright stars in our galaxies. 
When this structural data is combined with information from the coming generation of 
astrometric satellites (FAME, SIM and GAIA), we will finally have a dynamically detailed 
picture of our home, the Milky Way. 

Thanks are extended to Janet Weiland for assistance in installing the UIDL COBE data 
analysis package. The COBE datasets were developed by the NASA Goddard Space Flight 
Center under the guidance of the COBE Science Working Group and were provided by the 
NSSDC. 
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Fig. 1. — The local error in the J passband, with respect to galactic latitude, evaluated at 
each pixel location within a 7x7 pixel window. 

Fig. 2. — Map of the ratio of the smoothed to unsmoothed data (top), and the estimated 
local error <t 2 4o (bottom). The maps are Mollweide galactic projections for |6| < 30deg, on 
logarithmic scales, indicated by the grey scale bars. Blanked pixels are regions excluded 
from the present analysis (see text). Positive galactic longitudes are to the left of center. 

Fig. 3. — The relative emissivity of the dust disk as a function of temperature galactocentric 
radius. 

Fig. 4. — 240/im skymaps of the smoothed DIRBE data (top), the modeled emission (mid- 
dle), and the relative difference between the modeled and observed emission (bottom). The 
top two maps are on a logarithmic scale, while the bottom map is on a linear scale. 

Fig. 5. — The 240/zm emission profile for the data (X's) and model (diamonds) within 
0.17deg of the GP (6 = 0), on a logarithmic scale (log.D). 

Fig. 6. — 240/im emission profiles within 0.17deg of the indicated galactic latitudes. Scale 
and symbols are as in previous figure. 

Fig. 7. — The logarithm of the absolute FIR deviations, between the data and the model, 
as a function of galactic latitude. 

Fig. 8. — Skymaps of the observed J band emission (top), the predicted J band emission from 
the stellar distribution model (middle), and the relative difference map between model and 
data (bottom). Black pixels in the upper two maps show pixels removed from the data set, 
due to blanking out selected regions, such as the Galactic center, or rejected point sources. 
Upper two maps are on a log scale, while the bottom is on a linear scale. 

Fig. 9. — Same as previous figure, but for the K band. 

Fig. 10. — Observed (X's) and modeled (solid line) J and K band emission within 0.2 deg of 
the GP for the sheared spiral arms model. Vertical scale is in magnitude units (—2.5 log_D b )- 

Fig. 11. — J band emission as observed and modeled within 0.2 deg of the indicated latitudes. 
Scale, symbols and units are as in previous figure. 

Fig. 12. — K band emission as observed and modeled within 0.2 deg of the indicated latitudes. 
Scale, symbols and units are as in previous figure. 
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Fig. 13. — Surface density map of the dust, as inferred from the dust density model. Small 
black dot (upper center) shows the position of the Sun, which nearly lies on a small local 
feature, known as the Orion arm. Arms are incomplete on the side opposite the Sun due to 
incomplete HII data. 

Fig. 14. — Surface brightness map of the Milky Way in the K band. Bright dot (upper 
center) indicates the position of the Sun. 

Fig. 15. — Scale heights as a function of radius for the stellar disk (dashed line), the dust 
disk (solid line), and spiral arms (dot dash line), assuming R e = 8kpc. 

Fig. 16. — Arm - interarm density ratio for the stellar (dashed curve) and the dust spiral 
arms (solid curve). 

Fig. 17. — Predicted J band GP emission profiles at \l\ < 90deg for an axisymmetric model 
and three spiral models, m = 2 and m = 4 logarithmic spiral models, and a sheared HII 
spiral model, plotted against the data (X's). 

Fig. 18. — Same as previous, but for the K band. 

Fig. 19. — Schematic of the Galaxy showing the four spiral arms as mapped by HII regions 
and the dust (bold lines), the sheared arms in the K band (stars), and the arms in the 
two-arm logarythmic model for J and K band fit (dashed) and the K band fit alone (solid) 
(Drimmel and Spergel 2001). The HII spirals are incomplete on the opposite side of the 
Galaxy due to lack of data. 
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Table 1. Dust parameters. 



Parameter Symbol Value Uncertainty 



Disk component: 








density 


Po 


1098 MJy/sr/kpc 


41 


scale length 


h r 


2.26 kpc 


0.16 


base scale height 


ho,d 


134.4 pc 


8.5 


linear flair coefficient 


hi,d 


14.8 pc/kpc 


9.9 


flair radius^ 


n 


4.40 kpc 




240^m offset 


<?240 


1.07 MJy/sr 


0.15 


Local Arm: 








Radius at <f> = 


Rs 


8.001 kpc 


0.001 


pitch angle 


Ps 


7.33 deg 


0.45 


gap parameters: 


&t 


1.500 deg 






02/^1 


-0.643 


0.017 


half width 


W a /<j>i 


16.28 pc/deg 


0.70 


height 


z s /4>i 


11.9 pc/deg 


2.9 


central density 


Ps 


156 MJy/sr/kpc 


29 


emissivities 


k+/k- 


3.09 


0.52 



Spiral arms: 



emissivity' 1 ' 


kn 


1.00 




density 


Pa 


162 MJy/sr/kpc 


28 


cut-off radius 


r m 


6.71 kpc 


0.48 


cut-off scale length 


r a 


2.07 kpc 


0.22 


arm width coefficient 


Ca 


64.1 pc/kpc 


8.8 


initial scale height 1 * 




80.0 pc 




flair radius 


n> 


5.48 kpc 


0.13 


flair coefficient 


hi. a 


14.6 pc/kpc 2 


1.4 


Sag-Car reduction factor 


fr 


0.407 


0.036 


Warp: 








Radius that warp starts 




6.993 kpc 


0.046 


Phase of the warp 




-0.12 deg 


0.45 


amplitude coefficient 




72.8 pc/kpc 2 


7.0 



t Fixed parameters. 
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Table 2. Stellar parameters. 



Parameter Symbol Value Uncertainty 



Disk: 








normalization in J 




14.7 MJy/sr/kpc 


2.3 (4.2) 


normalization in K 


< 


11.6 MJy/sr/kpc 


1.7 (3.2) 


radial scale length 




2.264 kpc 


0.083 (0.19) 


scale height 


h* 


282.2 pc 


7.9 (20.) 


cutoff radius 


r c 


10.52 kpc 


0.34 



Spiral arms: 



amplitude in J 


Bj 


0.86 


0.25 


amplitude in K 


B K 


1.28 


0.24 


mean age in J 


TJ 


5.5 Myr 


3.6 


mean age in K 


TK 


17.6 Myr 


3.8 


arm width coeff. 


r* 


142. pc/kpc 


43. 


arm width coeff. 


dk 


143. pc/kpc 


69. 


corotation radius 


Rc 


6.66 kpc 


0.39 



Dust: 



V opacity 




0.0180 (MJy/sr)- 1 


0.0029 


emissivities 


K 


3.98 


0.35 (0.45) 






1.29 


0.10 (0.26) 






2.07 


0.37 (2.9) 


densities 


Ps 


121. MJy/sr/kpc 


23. (32.) 




Pa 


61. MJy/sr/kpc 


16. (26.) 


Miscellaneous: 








height of Sun 


Z Q 


14.6 pc 


2.3 


warp coefficient 


d w 


27.4 pc/kpc 2 


2.5 


J offset 


Qj 


-0.0684 MJy/sr 


0.0082 


K offset 


Qk 


-0.0744 MJy/sr 


0.0069 



Table 3. Dust models with variable constants. 







^0,a = 




ff, a = 




<t>i = 




Zq = 




VT = 




Parameter 


Units 


60 pc, 


100 pc 


O.5i? , 


O.6R 


1.0 deg, 


3.0 deg 


10 pc, 


20 pc 


-5.5K/R , 


-7.5K/i? Q 


multicolumnllDisk: 
























PO 


MJy/sr/kpc 


1122 


1112 


1057 


1096 


1120 


1172 


1018 


1060 


1084 


1067 


h r 


kpc 


2.21 


2.28 


2.26 


2.31 


2.19 


1.93 


2.52 


1.99 


2.22 


2.32 




pc 


139.0 


131.5 


132.4 


135.9 


135.3 


157.3 


135.0 


153.1 


136.7 


134.5 


n l,d 


J-"- / rV|JL 


16.5 


14.7 


15.3 


12.3 


20.9 


44.1 


0.8 


46.4 


13.5 


14.5 


Q240 


1\/T Tir / cr 

ivi jy / si 


1 04 


1 is 

1.1U 


1 07 


1 1Q 
1 . IO 


1 on 


n Q7 


1 42 






1 HQ 


r - * — i ttt i — a 

multicolumnllljocal Arm: 
























Ps 


deg 


7 QC 
/ .OO 


7 Qft 


ft fiQ 
0.00 


7 QQ 
i .00 


7 Qft 


ft Q7 


7 7/1 

/ . ( 4 


7 fiC 
/ .00 


7 17 


7 01 






-U.04o 


-U.OoD 


-U.0O4 


n ftftn 
-U.OOU 


-U.D04 


n ft/in 
-U.04U 


n ftftft 


n din 

-u.oiy 


-U.OoU 


n fte;7 
-U.Oo 1 


w B /4>i 


pc/deg 


16.43 


16.25 


16.65 


16.21 


16.88 


18.46 


15.85 


17.08 


17.09 


16.50 


7, /rh-i 


i"\ /■> / /loir 

pc/aeg 


12.0 


12.1 


11.8 


11.9 


17.5 


5.7 


9.5 


14.7 


11.8 


11.8 




MJy/sr/kpc 


161 


158 


160 


155 


223 


101 


152 


183 


152 


158 


h , Ih 

K,-\- j fi. — 




2 84 


2 82 


9 71 


3 1 ^ 




1 A7 


0.01 


9 HQ 


9 7fi 


z.ou 


multicolumnllSpiral arms: 
























Pa 


MJy/sr/kpc 


174 


130 


165 


160 


148 


87 


139 


104 


164 


162 


r m 


kpc 


6.36 


7.45 


6.56 


6.77 


7.16 


7.85 


7.19 


7.53 


6.55 


6.73 


ra. 


kpc 


2.23 


1.68 


2.15 


2.02 


1.81 


1.59 


1.78 


1.74 


2.14 


2.05 


Ca 


pc/kpc 


69.6 


57.4 


66.3 


60.7 


67.6 


81.4 


53.1 


81.5 


64.9 


64.1 


rt,e. 


kpc 


5.38 


5.71 


5.45 


5.46 


5.69 


5.38 


5.32 


5.34 


5.42 


5.48 




pc/kpc 2 


14.7 


14.5 


14.2 


14.9 


14.6 


11.3 


15.3 


11.0 


14.2 


14.5 


fr 




0.432 


0.384 


0.406 


0.404 


0.409 


0.518 


0.420 


0.440 


0.404 


0.406 


multicolumnllWarp: 
























-Rw 


kpc 


7.016 


6.973 


6.990 


7.000 


7.017 


7.053 


7.132 


7.041 


6.966 


6.996 




dcg 


0.18 


0.05 


0.02 


-0.30 


0.43 


1.28 


0.44 


0.57 


-0.14 


-0.11 


a„ 


pc/kpc 2 


75.9 


73.6 


72.7 


73.8 


77.5 


86.7 


91.8 


84.4 


70.3 


73.4 



Table 4. Alternative Stellar models. 







fro, a = 




ff,a = 




<t>i = 




Zq = 




VT = 




Parameter 


units 


60 pc, 


100 pc 


0.5-Rq, 


0.6R© 


1.0 deg, 


3.0 deg 


10 pc, 


20 pc 


-5.5K/R©, 


-7.5K/_R Q 


niulticolurririlllJisk: 


























ivi jy/ sr/ Kpc 


1O.0 


7 £. 




ZU.O 


O.D 


10. D 




ID. O 




Q A 




MJy/sr/kpc 


14.7 


6.2 


11.5 


15.8 


6.9 


12.4 


11.5 


13.1 


11.6 


7.5 


T* 


kpc 


2.063 


2.596 


2.264 


2.101 


2.564 


2.220 


2.264 


2.186 


2.264 


2.514 




pc 


290.5 


329.5 


283.3 


265.1 


318.0 


283.6 


279.3 


284.7 


283.2 


309.4 




kpc 


11.48 


10.66 


10.55 


10.45 


10.58 


10.17 


10.36 


10.51 


10.59 


10.56 


1 — • 1 . " 

multicolumnllbpiral arms: 
























Rt 




l.OU 


u.oy 


u.oo 


U. ( u 




U.tJU 


u.oo 


U.JO 


u.oz 


1 m 

l.Ul 


B K 




1.24 


1.43 


1.30 


1.07 


1.62 


1.04 


1.15 


1.76 


1.29 


1.49 


r J 


Mvr 


12.7 


11.4 


9.0 


7.1 


7.4 


16.3 


12.5 


9.2 


3.7 


7.3 


r X 


l\/f vr 
iviyr 


24.5 


13.7 


17.2 


23.2 


13.7 


20.4 


19.0 


17.9 


14.9 


15.4 


C J 


v"\ /" • / r\ e • 
J-"- / JV|JL 


158. 


198. 


135. 


143. 


145. 


79. 


208. 


72. 


132. 


158. 


C K 


pc/kpc 


O-LU. 


iyo. 


ioy. 




lOU. 


1 nn 


zuu. 




1 1A 


1 fi^i 

1UU. 




kpc 


5.67 


6.81 


6.69 


6.12 


6.72 


6.44 


6.33 


6.86 


6.95 


6.62 


multicolumnllDust: 




























4.16 


3.46 


3.64 


4.36 


3.53 


2.64 


3.94 


4.03 


3.62 


3.42 






1.46 


1.23 


1.34 


1.38 


1.22 


1.80 


1.19 


1.93 


1.31 


1.19 


^a 




1.77 


2.46 


2.80 


2.37 


2.94 


10.65 


2.65 


7.99 


2.84 


2.56 


Ps 


MJy/sr/kpc 


110. 


128. 


119. 


112. 


183. 


56. 


128. 


95. 


116. 


132. 


Pa 


MJy/sr/kpc 


98. 


53. 


59. 


68. 


50. 


8. 


53. 


13. 


58. 


63. 


multicolumnllMisc. : 


























(MJy/sr)- 1 


0.0177 


0.0149 


0.0182 


0.0194 


0.0152 


0.0235 


0.0170 


0.0225 


0.0182 


0.0155 


Qj 


MJy/sr 


-0.0441 


-0.0704 


-0.0682 


-0.0680 


-0.0691 


-0.0672 


-0.0715 


-0.0608 


-0.0673 


-0.0711 


Qk 


MJy/sr 


-0.0565 


-0.0781 


-0.0743 


-0.0715 


-0.0791 


-0.0737 


-0.0762 


-0.0669 


-0.0738 


-0.0790 


Z & 


pc 


14.9 


14.4 


14.8 


14.9 


14.8 


14.8 


10.0 


19.8 


14.8 


14.5 




pc/kpc 2 


25.0 


25.4 


26.9 


29.8 


26.1 


27.6 


33.3 


26.5 


26.2 


26.5 
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Table 5. Alternative Stellar models. 



Parameter units stellar flair r c = oo Qj,k = 0. k a = 1 /i* = 2/i , a 



Disk: 





MJy/sr/kpc 


12.2 


13.2 


10.2 


14.1 


22.8 




MJy/sr/kpc 


9.7 


10.3 


7.4 


11.1 


16.9 


r* 


kpc 


2.319 


2.255 


2.302 


2.273 


2.06 


h* 


pc 


292.3 


309.0 


306.30 


285.6 


253.2 


r c 


kpc 


11.10 


oo 


12.41 


10.54 


10.90 



Spiral arms: 



Bj 




1.06 


0.97 


0.69 


1.33 


0.14 


B K 




1.62 


1.24 


1.40 


1.64 


0.44 


Tj 


Myr 


4.0 


-13.2 


17.5 


7.5 


15.6 


TK 


Myr 


11.8 


13.3 


27.6 


14.3 


26.0 


r* 


pc/kpc 


155. 


266. 


136. 


173. 


95. 


r* 


pc/kpc 


151. 


262. 


413. 


164. 


328. 


Rc 


kpc 


6.34 


6.95 


5.96 


6.85 


6.44 



Dust: 



Ps 
Pa 



(MJy/sr)- 1 



MJy/sr/kpc 
MJy/sr/kpc 



0.0162 


0.0178 


0.0116 


0.0180 


0.0188 


3.95 


4.17 


4.06 


4.02 


3.38 


1.28 


1.35 


1.31 


1.30 


1.09 


2.17 


2.08 


2.06 


1.00 


1.59 


122. 


115. 


118. 


119. 


142. 


75. 


78. 


79. 


162. 


102. 



Miscellaneous: 

Z Q pc 14.8 13.2 14.8 14.7 14.5 

a w pc/kpc 2 24.5 19.4 22.0 29.1 25.6 

Qj MJy/sr -0.0731 -0.0585 0.0000 -0.0687 -0.0621 

Q K MJy/sr -0.0798 -0.0648 0.0000 -0.0740 -0.0620 
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Table 6. Alternative Spiral models. 



Parameter 


units 


axisymmetric 


m=2 log 


m=4 log 


sheared 


Disk: 












n 

m 


MJy/sr/kpc 


14.9 


15.1 


11.8 


14.7 


n 

Vk 


MJy/sr/kpc 


12.6 


12.5 


9.9 


11.6 


r* 


kpc 


2.257 


2.255 


2.386 


2.264 


ft* 


pc 


266.8 


274.1 


286.0 


282.2 


r c 


kpc 


10.70 


10.90 


10.03 


10.52 


Spiral arms: 












sheared: 






















0.86 


B K 










1.28 


tj 


Myr 








5.5 


TK 


Myr 








17.6 




pc/kpc 








142. 




pc/kpc 








143. 


Rc 


kpc 


- 


- 


- 


6.66 


logarithmic: 












B 






1.38 


0.42 




P 


deg 




25.4 


15.8 




i?o 


kpc 




5.38 


2.56 




< 


pc/kpc 




138. 


99. 




Dust: 












k+ 
n/ + 




2.65 


3.61 


3.53 


3.98 


fc_ 




0.86 


1.17 


1.14 


1.29 






1.00 


2.19 


11.25 


2.64 


Ps 


MJy/sr/kpc 


181. 


133. 


136. 


121. 


Pa 


MJy/sr/kpc 


162. 


74. 


14. 


61. 


Miscellaneous: 












Ky 


(MJy/sr)- 1 


0.0143 


0.0172 


0.0158 


0.0180 


Qj 


MJy/sr 


-0.0533 


-0.0633 


-0.0640 


-0.0684 


Qk 


MJy/sr 


-0.0744 


-0.0760 


-0.0814 


-0.0744 


Z G 


pc 


14.7 


14.8 


14.9 


14.6 


a w 


pc/kpc 2 


25.3 


24.0 


28.4 


27.4 


X 2 for \b\ < 3deg 




65502 


51547 


53047 


51452 



-41 - 



-42- 




-20 30 

latitude 



-43- 




-45 - 




-46- 



3.5 




-47- 




1DD D -1DD 1DD D -1DD 

gal long gal long 



-48- 



T 




J I I I I I I I I I l_ 

-20 20 

latitude 



- 50 - 






- 52 - 




Lj i i I i i i i I i i i i I i i i -L i" i i I i i i i I i i i i I i i 

1DD D -1DD 1DD D -1DD 



gal long gal long 



- 53 - 




-54- 




- 58 - 



i i i | i i i i | i i i i | i i i i i i | i i i i | i i i i | i i r 




j i i I i i i i I i i i i I i i L 



50 D -50 50 D -50 

gal long gal long 



- 59 - 
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gal long gal long 



